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Abstract 

One of the most promising possibilities may be the appearance of quark mat- 
ter in astrophysical phenomena in the light of recent progress in observations. 
The mechanism of deconfmement is not well understood, but the thermodynam- 
ical aspects of the hadron-quark (HQ) phase transition have been extensively 
studied in recent years. Then the mixed phase of hadron and quark matter 
becomes important; the proper treatment is needed to describe the HQ phase 
transition and derive the equation of state (EOS) for the HQ matter, based on 
the Gibbs conditions for phase equilibrium. We here adopt a EOS based on 
the baryon-baryon interactions including hyperons for the hadron phase, while 
we use rather simple EOS within the MIT bag model for the quark phase. 



*E-mail address: nobutoshi.yasutake@it-chiba.ac.jp 



2 



Nobutoshi Yasutake,Tsuneo Noda et al. 



One of the interesting consequences may be the appearance of the inhomoge- 
neous structures called "pasta", which are brought about by the surface and 
the Coulomb interaction effects. We present here a comprehensive review of our 
recent works about the HQ phase transition in various astrophysical situations: 
cold catalyzed matter, hot matter and neutrino-trapped matter. We show how 
the pasta structure becomes unstable by the charge screening of the Coulomb 
interaction, thermal effect or the neutrino trapping effect. Such inhomogeneous 
structure may affect astrophysical phenomena through its elasticity or thermal 
properties. Here we also discuss some implications on supernova explosion, 
gravitational wave and cooling of compact stars. 

1. Introduction 

The mechanism of the hadron-quark (HQ) deconfinement is one of hot topics in 
modern physics. Theoretically, numerical simulations based on the lattice gauge 
theory have been done to estimate the critical temperature and figure out proper- 
ties of the phase transition pQ. However, at present, these attempts do not provide 
definite results for the deconfinement transition in high density region. On the 
other hand, there have been done many theoretical studies by using the effective 
models of QCD[2l [3]. Therefore, one of the most hopeful strategies to get insight 
into this mechanism is to find or predict some signals about QCD phase transition 
in the astrophysical phenomena such as supernovae, black hole formations, mergers 
of neutron stars, cooling of compact stars, etc., in the light of recent theoretical 
and observational progress. Then the equation of state (EOS) is a basic element 
to discuss the implications on above astrophysical phenomena. There have been 
proposed and discussed various types of scenario concerned with the HQ deconfine- 
ment transitions at high-density region, but it is still unclear, even, whether the 
phase transition is the cross over or the first order [3]. Only assumption in this 
review is that the phase transition is of the first order as suggested by many model 
studies [2]. One of the direct consequences of this assumption is the emergence of 
the hadron-quark (HQ) mixed phase during the phase transition. In this review we 
first discuss the thermodynamical aspects of the HQ transition, and present some 
implications on astrophysical phenomena subsequently. 

Historically, the mixed phases have been often treated very naively by applying 
the Maxwell construction to get the EOS in thermodynamic equilibrium. In 1990 
Glendenning has pointed out that the phase equilibrium in multi-component system 
or in the case of plural chemical potential corresponding to the conserved quantities 
such as in neutron star matter must be considered based on the more general Gibbs 
conditions [51 [6] . He emphasized that the neutron-star matter should be treated as 
abinary system specified by two conserved quantities, electromagnetic charge and 
baryon numbers. Thus the pressure becomes no more constant as in the Maxwell 
construction, but still increases in the mixed phase. He also demonstrated that the 
mixed phase spreads in a rather wide density region within a bulk calculation, where 
two kinds of semi-infinite matter are considered in phase equilibrium, and the pres- 
sure balance, chemical equilibrium and the global charge-neutrality conditions are 
imposed, discarding the Coulomb interaction. It is important to see that particle 
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fractions in each phase are no more the same for given average particle-fractions. 
Generally any phase transition in a system with two or more conserved charges must 
be non- congruent^ . However, the Coulomb interaction is discarded in this calcula- 
tion and we can by no means estimate how important it is, or the surface energy can 
not be taken into account a priori. If we consider a more realistic situation, we must 
consider inhomogeneous matter with various geometric structures called " pasta" , 
taking into account the finite-size effects such as the Coulomb interaction and the 
surface tension. The pasta resembles the nuclear pasta in liquid-gas transition at 
subnuclear densities [H El [91 [TUl [HI Q2]. Heiselberg et al. have demonstrated that 
such finite-size effects are important for the mixed phase in the context of the HQ 
transition [13] . Treating the pasta structure as in ref. [13], they simply considered 
uniform particle densities in both phases. They showed a possibility of that the 
finite-size effects may largely restrict the region of the mixed phase. Subsequently, 
Voskresensky et al. have shown that the rearrangement of particle densities and the 
screening effect for the Coulomb interaction are important, and the large surface 
tension gives the mechanical instability of the geometrical structures in the pasta 
phase [15] . Using the effective model for hadron matter and the MIT bag model for 
quark matter, Endo et al. explicitly demonstrated these features [16] , Following this 
idea, we have studied the pasta structures brought by the kaon condensation and 
the HQ deconfinement transition \17\ [T6| [T8] . More recently we have extended our 
framework to deal with the finite temperature case or the neutrino-trapping matter, 
which may be relevant to supernova explosions or neutron star mergers [19\ 120] , 

In this review we first present thermodynamic concepts of the non-congruent 
transition of the HQ phase transition (Section [2j. In the section we will also present 
some results about the phase transition, where the effects of finite temperature or 
neutrino trapping are also discussed. In Section [3]] is discussed some astrophysical 
implications of the HQ phase transition. Section [T] is devoted to summary and 
concluding remarks. 

2. Phase Equilibrium in Thermodynamics on Hadron- 
Quark Phase Transition 

In this section we discuss the emergence of the pasta structure associated with the 
HQ deconfinement transitional For our purpose we need both EOS's of hadron and 
quark matter as realistically as possible. As we mentioned in the last section, no 
one knows how to exactly calculate the HQ phase transition at high density region. 
The studies by using the effective models of QCD such as the MIT bag model or 
the NJL model have been actively done instead. We, hereafter, use the MIT bag 
model for simplicity. 

We adopt a realistic EOS for the hadron matter based on the microscopic 
nucleon-nucleon (NN), and nucleon-hyperon (NY) interactions. Although there 

1 When one dare to apply the Maxwell construction in such systems, one must violate the Gibbs 
conditions by ,e.g., imposing local charge neutrality. 

2 We, hereafter, dismiss the possibility of color superconductivity [21] . 
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are many models to describe the EOS with hyperons, they have many uncertainties 
which come from the uncertainty of NY and YY interactions. However, nowadays, 
the potentials of the baryon-baryon interactions have been directly calculated by 
the lattice QCD calculations [22]. We hope that they provide realistic NY and YY 
interactions in near future. Motivated by such background, it becomes a very hot 
topic to calculate the EOS from the baryon-baryon interaction directly. Here, we 
show the EOS by the theoretical framework of the nonrelativistic Brueckner-Hartree- 
Fock (BHF) approach [23] based on the microscopic NN and NY interactions. The 
BHF calculation is a reliable and well-controlled theoretical approach for the study 
of dense baryon matter. 

- Hadron Phase 

" BHF theory " 
EOS based on NN, and NY interactions. 



1 1 

- Mixed Phase 

" pasta structures " 

• chemical equilibrium 

• baryon/fepton number conservation 

• charge neutrality 

• surface energy/charge screening 

-Quark Phase 

"Bag model, NJL model, etc." 
Effective theory. 



Figure 1. Schematic picture of our strategy on the study about the HQ phase 
transition for the astrophysical implications. 

In Fig. [H we sum up our strategy to construct EOS including the HQ phase for 
the astrophysical implications. Since our goal is to figure out the realistic properties 
of the HQ phase transition, we chose one of the most realistic ways to know the HQ 
mixed phase at high-density region; NY, YY interactions will be clear in the future 
for the hadron phase. 

Then we consider the phase-equilibrium by including the finite-size effects. The 
Coulomb interaction and the surface energy are taking into account in a consistent 
way: charge screening effect is automatically took into account by solving the Pois- 



We assume 

"the first order phase transition" 
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son equation, but the surface energy at the interface of hadron and quark matter is 
treated in a phenomenological way by using the surface tension parameter. Thus, 
the largest ambiguity in our study is the uncertainty in the quark models. Hence, 
one of the main purposes of our study is to avoid this ambiguity as much as possi- 
ble by comparing our results with astrophysical phenomena. Here we sum up our 
method, and report our current status. 

2.1. Hadron Phase; the Brueckner-Hartree-Fock theory 

BHF theory is well-known to describe hadron matter using the realistic baryon- 
baryon interaction. In this review we do not explain the details of the procedure, but 
they can be found in Refs. [231 [Ml US [26] . We adopt the Argonne Vis potential [27] 
for NN interaction, and semi-phenomenological Urbana UIX nucleonic three-body 
forces [29] and the Nijmegen soft-core NSC89 NY potentials [30]. Unfortunately, 
there are not reliable potentials for the hyperon-hyperon interaction now. Therefore, 
we neglect them here. 

To take into account the effects of finite temperature, we adopt Frozen Cor- 
relations Approximation [31, 32, 33J. The correlations at finite temperature are 
assumed to be the same with the ones at zero temperature. It is found to be a good 
approximation at finite temperature as Refs. [3T | [32 ]. [55] . 

In the following, we briefly describe our framework. First, we get the chemical 
potential m from the number density Hi, 

9 f°° 



(2vr) 2 



o 



1 

h{p) = exp{ta-^)/T} + l' (2) 

where £j and fi{p) are the single-particle energy and the Fermi-Dirac distri- 
bution function, respectively, whereas subscript i denotes the particle species, 
i = n,p, A, £~. We set each degeneracy factor g = 2, and adopt each mass as 
m n = m p = 939 MeV, m A = 1115.7 MeV, and m E - = 1197.4 MeV. Note that e { 
includes the interaction energy {7, as well as the kinetic energy |31[ 134] . 



/ m 2 +p 2 + U ._ (3) 

Finally, we get the free-energy density T as follows, 
a r°° i I 



roc i i 

J y m| + p 2 fi(p) 4irp 2 dp + - UiUi 



(2tt) 3 . 

where sh is the entropy density calculated from 



Ts H , (4) 



SH = - E 7^)3 J Q {Mp) ln &(p) + C 1 - h(p))H^ - Up))} Wdp. (5) 



The total pressure for the uniform hadron phase is given by 

E 



Ph = y~] mm - Th- (6) 
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2.2. Quark Phase; Bag model 

Next, we explain how to treat the quark phase based on the thermodynamic bag 
model. In this model, we get the number density uq, the pressure Pq, and the 
energy density eg, 



E (2^)3 ^Jm 2 q +P 2 f q {p)Wdp + B, (8) 

S Q = -Et^W^ {/ g (p)ln/ 9 (p) + (l-/ g (p))ln(l-/ g (p))}47rp> (9) 
Jq = e Q -T SQ (10) 

for uniform quark matter, where f q (p) is the Fermi-Dirac distribution function of 
the quark q (= u,d,s), m q its current mass, and B the energy-density difference 
between the perturbative vacuum and the true vacuum, i.e., the bag constant is set 
to be 100 MeV-fm -3 . The baryon density n# is given by (n u + rid + n s )/3. 

Since there are still many uncertainties for the vacuum property, we here also 
adopt another type of B by the density dependent bag model [33] • In this model, 
we set the vacuum energy B(ng) as 



B(n B ) = -Boo + (-Bo - -Boo) exp 



n 



with B^ = 50 MeV-fm~ 3 , B = 400 MeV-fm" 3 , and (3 = 0.17. Here, we set the 
saturation density as rig = 0.17 fm -3 . This model has been proposed in Ref. [33], 
on the basis of experimental results obtained at CERN SPS on formation of a 
quark-gluon plasma. 



Here, we show free energy per baryon F/A as a function of baryon density for 
hadron matter and quark matter at different temperature (T = 0, 10, 30, 50 MeV) 
in Fig. [2j We adopted the thermodynamic bag model (B = 90 MeV/fm 3 ) as the 
quark matter for left panels, but the density dependent bag model for the right 
panels. The lower panels show the case with neutrino trapping, Yj, = 0.4, while 
upper panels for without. Note that the values of F/A of hadrons are lower than 
the ones of quark matter at neutrino trapped case at T = 0, 10, 30 for any density. 
Therefore it is not appropriate for our assumption that the HQ transition of the 
first order occurs at high density region. Hence we adopt the density-dependent bag 
model for the lepton trapping case in this article. 

We note that the bag model should be quite simple to describe the quark matter. 
We will update the quark model to more sophisticated models such as NJL model [35J 
EE1EZ], pNJL model [381 EH 00 SD S21 S3] , or Dyson-Schwinger model [H], etc.. 
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Figure 2. Figure from Ref . [33] . Free energy per baryon as a function of baryon 
density of beta-stable hadron matter (solid lines) and quark matter (dashed lines) 
with a bag constant B = 90 MeV/fm 3 (left panels) or a density-dependent bag 
parameter (right panels) with (lower panels) and without (upper panels) neutrino 
trapping at different temperatures T = 0, 10, 30, 50 MeV (upper to lower curves). 

2.3. The Maxwell Construction and the Gibbs Conditions 

2.3.1. Hadron-quark Phase Transition as a Pure System Under the 
Maxwell Construction 

In this subsection, we introduce the mixed phase to discuss the phase equilibrium 
between hadron and quark phases. We can expect the HQ transition to appear in 
various circumstance, supernovae, protoneutron stars, cold neutron stars or neutron- 
star mergers, where temperature may become high or neutrinos may be trapped. 
So we consider here a general situation, where exist hadrons such as nucleons and 
hyperons and leptons such as electrons and neutrinos. 

First, we consider the HQ phase transition without considering the finite-size 
effects which is called the bulk calculation. The naive treatment of the phase equi- 
librium in the first-order phase transition should be the application of the Maxwell 
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construction. In this treatment for neutron-star matter we implicitly assume that 
electric charge is always locally neutralized in each phase and lepton number is equal 
to each other; only baryon number is different between two phases. So the Maxwell 
construction means that the matter is regarded as a pure system, discarding the 
electromagnetic charge and lepton number. 

Relaxing the local charge-neutrality condition and equality of lepton number, we 
must impose the Gibbs conditions in the proper treatment of the HQ transition in 
neutron stars. In this treatment, we impose the balance of the chemical potentials of 
quarks and baryons as well as the pressure balance and thermodynamical equilibrium 
between two phases. The latter two conditions are obvious, 

P H = P Q , T h = Tq, (12) 

where Pjj (Pq), Th (Tq) are pressure and temperature of hadron (quark) phase. 
The former conditions are reduced to only one within the Maxwell construction for 
the baryon- number chemical potential hb H , 

/xe(hadron phase) = //b (quark phase), (13) 

where the baryon chemical potential in the hadron (quark) phase is given by hb = 
Hn (hb = Hu+%Hd)- Here, Hm Hm and Hd is neutron, u-quark, and d-quark chemical 
potential, respectively. 

We also assume the beta equilibrium among the particles as 

1 2 Q 1 1 Q 

Hu = gMB + ^H C ' y-d = Ms = gMB - gMc' ( 14 ) 

for the quark phase, 

Hn = A* A = HB, fJ, p = fJ,B + He i M£~ + Mp = 2 MB> (15) 

for the hadron phase, where the quark chemical potentials are shown as Hu,d,s f° r 
u, d, s-quarks, the hadron chemical potentials HnpA,B- f° r n,p,A,T,~. In these 
equations, we have also introduced the two kinds of variables to understand the 
mechanism of the HQ phase transition systematically; the lepton chemical potential 
fi^ and the charge chemical potential fJ-Q are given as 

£} Q) =^ Q \ M? (Q) =M? (Q) -Mf (Q) , (16) 

in hadron (quark) matter. Here the lepton chemical potentials are shown as He 
and H?} f° r electrons and electron type neutrinos, respectively. Note that hl 
and He are not equal to each other between two phases; Hl ^ Ml an d He ^ Mc- 
Under this construction, we assume only the equality of baryon chemical potentials 
Hb between two phases shown as Eq, (fT3|) . Hence the matter is regarded as a pure 
system in this treatment. As a result, the EOS has a density jump in general. 



3 Another version of the Maxwell construction may be possible [45], where the effects of the 
lepton chemical potentials on fiB are taken into account. 
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2.3.2. Hadron-quark Phase Transition as a Ternary System Under the 
Gibbs Conditions 

So far the charge equilibrium is completely discarded between two phases: the above 
assumption uses the charge neutrality condition in the local form, while it must be 
attained as a whole system. Consequently, the Maxwell construction arises as a 
consequence of further strong constraint. In such treatment, the lepton number is 
also assumed to be conserved locally. In the proper treatment of the HQ transition 
we must consider both phases as ternary systems specified by the three conserved 
quantities, electromagnetic charge, lepton number, as well as baryon number. 

In other words, we must add two chemical equilibrium-conditions into Eq.(|13p 

as 

HQ HQ HQ 

Here, there are three identical components; the baryon number (which corresponds 
to hb), charge (fJ>c), and lepton number The Maxwell construction is then 

meaningless and we must use the Gibbs conditions to describe the phase equilibrium 
in a thermodynamically consistent way. Hence, the EOS dose not have density jump 
although we assume that the phase transition is of the first order. 

Since we impose the global charge neutrality (GCN), the total charge as well 
as baryon-number are shared by the two phases, conservation should be satisfied as 
following; 

n B = Xni + (l-X)r>>B > ( 18 ) 
Y c n B = X ^nJ + (l-x)y c H nf > (19) 
Y L n B = X Y? n Q B + {l- X )Y? nf , (20) 

where nB,n B ,n B , and \ are respectively the total baryon number density, the 
baryon number density in hadron matter, the effective baryon number density in 
quark matter, and the volume fraction which indicates the volume of the quark 
phase divided by the total volume. Here, is one third of the total quark number 
density. The values of Yq and Yl are the total fraction of positive particles and 
leptons, whereas Y®,^ and Yq^ denote the fractions of positive particles (leptons) 

per baryon in quark matter, and in hadron matter. Note that Yeas = YM^ are 

assumed in the previous treatment (fl~3|) instead of /w^) = ^C(L) • 



2.3.3. Hadron-quark Phase Transition as a Binary System with the Lo- 
cal Charge Neutrality 

Recently the hybrid types of the Maxwell construction and the Gibbs conditions 
have been suggested, depending on the physical situation [35]. One of them is the 
Gibbs conditions assuming the local charge neutrality (LCN), = Yq = 0. In 
this treatment, the chemical equilibriums are described as follow; 

H Q H Q /m\ 

fJ-B = Mb> Ml = Mi- (21) 
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We shall see LCN is turned out to be achieved in some cases as in the protoneutron 
stars (subsection 12. 8, p . 

In the last of this subsection, we show some examples of EOS adopted the Gibbs 
conditions for the HQ mixed phase in FiglSJ where we compare both conditions of 
GCN and LCN. The left panel is at finite isentropic case s = 2 [/ks], and Yl = 0.35 
assuming the matter in protoneutron stars. One can see that the proper treatment 
of the Gibbs conditions (GCN) gives almost the same result as that in the hybrid 
treatment of the Gibbs conditions with a further constraint, LCN, in this case. 
The right panel is at zero temperature without neutrinos assuming the case in the 
neutron stars. Note that, without neutrinos, the Gibbs conditions with LCN give 
completely the same result with the Maxwell construction. Consequently, the EOS 
has a density jump as shown in the right panel in Figj3] while pressure smoothly 
increases in GCN. We also show the uniform EOSs of quark matter and hadron 
matter in these figures. 




Figure 3. Equation of sate considering the mixed phases imposed the Gibbs con- 
ditions assuming GCN and LCN, which is labeled or 'GCN' and 'LCN'. For com- 
parison, we also shows the uniform matter of quark matter (labeled as 'quark') and 
hadron matter. Here, the EOS labeled as 'hyperon' shows the hadron matter includ- 
ing hyperons, while the EOS labeled as 'nucleon' is without hyperons. Left panel 
shows the EOS at isentropic, s = 2[/ks], and at Yl = 0.35. Right panel is the same 
with left one, but for at zero temperature without neutrinos. 

2.4. Finite-size Effects and "Pasta" Structures 

We have presented general features of the non-congruent transition in a clear fashion 
with a help of the bulk calculation, discarding the Coulomb interaction and surface 
effects called the finite-size effects. However, the bulk calculation is insufficient to 
figure out the essential aspects of the non-congruent transition. Actually we can 
see the emergence of imhomogeneous structure with various geometrical shapes, 
called pasta (see Fig. H]), by the cooperative effect of the Coulomb interaction and 
the surface energy, which is similar to the nuclear pasta at subnuclear densities 
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El El [ini El 02] • The charge screening effect, which is a typical many-body effect 
in the presence of the Coulomb interaction, plays an important role in this context to 
make the pasta structures unstable, which is called mechanical instability [15]. The 
evaluation of the free energy of the pasta structures can be carried out by using the 
Wigner-Seitz approximation; dividing the whole space into equivalent unit cells, we 
solve the coupled equations for density profiles of particles and the Poisson equation 
in a selfconsistent way within the unit cell |46j . We further use the Thomas- Fermi 
(local density) approximation for each particle. 



uniform droplet rod slab tube bubble uniform 




Figure 4. Schematic picture of "pasta" structures. When red phase appears in the 
blue phase, the red phase form droplets in the sea of blue. With increase of the 
volume fraction of red phase in the blue phase, the shape of the red phase changes 
from droplet to rod, slab, tube, and to bubble before pure red phase appears. 

Assuming various geometrical structures, we must numerically search the mini- 
mum of the free energy 

F = E v + F e + F Ve + E c + E s (22) 

changing the volume ratio of quark matter to hadron matter, and the size of Wigner- 
Seitz cell Rw- Here the bulk energy Ey gives a main contribution, 

E V = I dr 3 F H [ni]+ [ dr 3 T Q [n a ], (23) 
JV H JVq 

where Th and Tq comes from the kinetic energy and strong interaction. The free 
energies of electrons, F e , and neutrinos, F Ue , are simply calculated by the use of 
the Fermi-Dirac distribution functions. Eq denotes the interaction energy due to 
the Coulomb force. The treatment of the Coulomb potential may deserve some 
discussions, since it includes a delicate problem and sometimes has given rise to a 
confusion |15j . First of all we must notice that charge chemical potential /Uc* itself 
is not physical quantity since it is gauge dependent in the presence of the Coulomb 
interaction Vcoul- Only their combination of nc~ ^Coui is gauge invariant. Consider, 
for example, electrons without any neutrino. Then the electron chemical potential 
H e is equal to fj,c, so that the Gibbs condition means [i e equals to each other in 
the mixed phase. On the other hand, the electron number density n e should be 
represented as (/U e — Vcoul) 3 / (3tt 2 ) in the gauge invariant way. Thus if Vcoul spatially 
varies in the mixed phase, each phase can have different number of electrons even if 
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He is the same in the both phases. In the extreme case Vcoul m &y resemble the step 
function and we have a simple picture of uniform electron density in both phases 
with different numbers. Note that there arises no such problem in the absence of the 
Coulomb potential as in the bulk calculation: the difference of the electron number 
directly reflects in the difference of the electron chemical potential in this case. 

Here, we do not consider all kinds of anti-particles for simplicity. After the 
minimization of the free energy, we can get the thermodynamic quantities, such as 
pressure, entropy, etc.. The notation Eg stands for the surface energy which comes 
from a sharp interface between the hadron and quark phases with a fixed surface 
tension a. Although the surface tension of the HQ interface is poorly known, some 
theoretical estimates based on the MIT bag model for strangelets [47] and lattice 
gauge simulations at finite temperature [48] have suggested a range of a rj 10- 
100 MeV • fm~ 2 . We will discuss the effects of the surface tension on the pasta 
structures in subsection 12.6.1 
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Figure 5. Density profiles and the Coulomb potential Vc within a 3D (quark 
droplet) Wigner-Seitz cell of the mixed phase at = 0.4 fm -3 . The cell radius 
and the droplet radius are Rw = 26.7 fm and R = 17.3 fm, respectively. 

Fig. [5] shows an example of the resulting density profiles in the droplet phase 
for tib = 0.4 fm -3 at zero temperature |18j . One can see the non-uniform density 
distributions of particle species together with the non-vanishing Coulomb potential. 
The quark phase is negatively charged, so that d and s quarks are repelled to the 
phase boundary, while u quarks gather at the center. The protons in the hadron 
phase are attracted by the negatively charged quark phase, while the electrons are 
repelled. 

2.5. Hyperon Suppression in the Mixed Phase 

Though we allow the freedom of hyperon mixture, hyperons do not actually come 
into the mixed phase. The suppression of hyperon mixture in the mixed phase is due 
to the fact that the hadron phase is positively charged. As shown in Fig.0 hyperons 
appear in charge-neutral hadronic matter at low density to reduce the Fermi 
energies of electron and neutron. Without the charge-neutrality condition, on the 
other hand, there is symmetric nuclear matter at lower density and hyperons will be 
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mixed there at higher density due to the large hyperon masses. Generally speaking 
hyperons are hardly mixed in the positively charged matter. Thus, the mixture 
of hyperons is suppressed in the mixed phase where the hadron phase is positively 
charged. 




0.0 0.5 1.0 0.0 0.5 1.0 

« B [fm" 3 ] « B [fm" 3 ] 

Figure 6. (a) Particle fractions of neutral matter with electrons, (b) The same 
quantity for charged matter without electrons, the low-density part of which corre- 
sponds to symmetric nuclear matter. 

2.6. Effect of the Surface Tension 

Let us consider the role of the surface tension on the mixed phase. As already 
mentioned, the surface tension in the interface of the hadron and quark phases is 
poorly known. If one uses a smaller surface tension parameter a, the energy gets 
lower and the density range of the mixed phase gets wider. The limit of a = MeV- 
fm~ 2 leads to a bulk application of the Gibbs conditions without the Coulomb and 
surface effects. On the other hand, using a larger value of a, we find that the scale 
of the inhomogeneity becomes large and the EOS gets closer to that of the Maxwell 
construction case. Beyond a limiting value of a ~ 65 MeV • fm -2 the structure of 
the mixed phase becomes mechanically unstable [15J: for a fixed volume fraction 
(R/R\y) 3 the optimal values of R and Rw diverge and local charge neutrality is 
recovered in the mixed phase, where the energy density equals to that of the Maxwell 
construction (see Fig. [7j). 

Note that these roles of surface tension on the EOS is remarkable for the matter 
without neutrinos, i.e. a binary system, where the pasta structure is rather stable. 
If we consider the neutrino contamination, the matter should be treated as a ternary 
system. As we will discuss at subsection 12.84 the EOS could not be close to the 
one of the Maxwell construction, but to the one with the Gibbs conditions (LCN) 
depending on the neutrino fraction. 

2.7. Thermal Effects 

In following we discuss the thermal and neutrino trapping effects on the HQ phase 
transition. In the left panel of Fig. [8] we show the dependence of the isothermal free- 
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Figure 7. Droplet radius (R) dependence of the energy per baryon for fixed baryon 
density ng = 0.35 fm -3 and different surface tensions. The temperature T is zero 
for all cases. The quark volume fraction (R/R\y) 3 is fixed for each curve. Dots on 
the curves show the local energy minima. The black line shows the energy of the 
Maxwell construction case. 

energy per baryon on the cell size in Eq. (|22p for several temperatures at tib = 2 no, 
using the constant bag constant and the surface tension as B = 100 MeV-fm -3 and 
a = 40 MeV-fm -2 , respectively [19] . Note that the free energy per baryon AF/A 
is normalized value by subtracting the one at infinite Rw, for the comparison of 
thermal effects. Here we assume the geometrical structure is the droplet and do not 
take into account neutrinos to see the genuine thermal effects. The quark volume 
fraction (R/Rw) 3 is fixed to be the optimal value at T = MeV for each curve. 
From the this panel, we can see that the minimum point is shifted to a larger value 
of R as temperature is increased, and eventually disappears for T > 60 MeV. The 
temperature dependence of the free energy comes from the Coulomb energy, the 
surface energy and the correlation energy. The correlation energy is given by the 
difference in the bulk energy Ey in Eq. (|22[) [15] , defined by 

£corr = E v (m(r)) - E v {n{' H ). (24) 

The rearrangement of the particle densities also induces a change in the volume 
energy Ey, which is composed of kinetic energy and strong-interaction energy. In 
the above calculation Ey itself is intrinsically large but independent on R, since we 
have assumed uniform densities Th' 11 in each phase. Once the Coulomb interaction 
is properly taken into account, particle densities are no more uniform to produce a 
R dependence in Ey. The detailed analysis showed that E cmr gives negative values 
and behaves i? -1 for large R. By comparing the temperature dependence of these 
contributions, we can see that the correlation energy is primarily responsible to the 
behavior of the minimum point. It means that the mixed phase becomes less stable 
as temperature is increased [19] . 

In the right panel of Fig. [8j we show the EOS (pressure vs baryon density) for 
the neutrino-free case at T = 30 MeV. Clearly, the pressure becomes to close to 
the one given by the Maxwell construction at finite temperature, since the mixed 
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phase becomes unstable and its density regime is thereby largely restricted. The 
mechanical instability of the pasta structure also arises in other geometrical shapes, 
which means that both R and R\y go to infinity with their ratio fixed. Thus we 
recover the picture of the phase equilibrium of two bulk matters, where the surface 
energy is irrelevant and the Coulomb energy vanishes due to the achievement of 
the local charge-neutrality. These features indicate that the EOS resembles the one 
given by the Maxwell construction. 
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Figure 8. Left panel shows that the radius R dependence of the free-energy per 
baryon for the droplet case at ub = 2 no and different temperatures. The free 
energy is normalized by the value at R — > oo. The filled circles on each curve shows 
the energy minimum. The results are for B = 100 MeV-fm -3 , a = 40 MeV-fm -2 . 
The right panel shows the EOS of the mixed phase (thick dots) in comparison with 
pure hadron and quark phases (thin curves) at T = 30 MeV. We also show, for 
comparison, the mixed phase by the Maxwell construction by thin solid line. 



2.8. Neutrino Trapping Effects 

In Fig. [9j we show the dependence of the free energy on neutrino trapping in the 
mixed phase [20] . Here, we adopt the density-dependent bag, and fixed the quark 
volume fraction (R/Rw) 3 at Y Ue = 0.01 for each curve. For Y Ve > 0.1, the minimum 
point disappears. The right panel shows the each contribution to the free energy. 
From this figure, we can see that both the correlation energy and the Coulomb 
energy mainly contribute to the behavior of the minimum point in the presence of 
neutrinos, while it is a only the correlation energy for the temperature dependence. 

To elucidate the neutrino-trapping effect in the mixed phase, we show the density 
profiles and the Coulomb potential for the slab case in Fig. [TUl For comparison, 
we use the same cell-size, and fix R\y here. High neutrino fraction enhances the 
number of electrons to satisfy the chemical equilibrium. As a result, the Coulomb 
potential drastically changes in the neutrino trapping case as shown in this figure. 

In the left panel of Fig. \TT\ we show the relation of the pressure and density 
for PNS matter (T = 30MeV and Yl = 0.4) as a result of above calculations. For 
comparison, we also check the phase transition under the Gibbs conditions with 
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Figure 9. Left panel is the same as the left panel of Fig. [8l but for different neutrino 
fractions Y Ue at a given temperature, T = 10 MeV. In this case, we set baryon density 
as hb = 2.5 no- The thick lines of the right panel are the case of Y Ue = 0.15, the thin 
lines Y Ue = 0.01. The right panel shows each component of the free energy (AF/A); 
e.g. the free energy, the Coulomb energy (Eq/A), the surface energy (Eg/ A), and 
the correlation energy (E COTV /A). Here they are shown per baryon. 



GCN and LCN. As shown in the left panel of Fig. [U there is no minimum point 
of the free energy for PNS matter, which means the crystalline structures of pasta 
are broken there to form gan amorphous stateh composed of two species of matter; 
the amorphous state should take a bicontinuous structure with a complex pattern, 
depending on the kinetics of its formation. In such case, the surface tension and the 
Coulomb interaction do not work to make any pasta structure (see the right panel 
of Fig. [9}, so that the mixed phase may resemble the amorphous phase of quarks 
and hadrons. For numerical reason, we can not calculate R\y = oo. Instead, we set 
Rw = 100 fm for amorphous state. 

From the left panel of Fig. [TTJ and the above discussion, the amorphous state of 
quarks and hadrons gives almost the same result with the bulk Gibbs calculation 
with LCN. To understand this reason, we also show the density profile and the 
Coulomb potential 3no assuming the rod structure for the same PNS matter in the 
right panel of Fig. [TTJ From this figure, we can see that the shape of the Coulomb 
potential looks like a step function. This is because the size of Rw is very large 
(Rw = 100 fm); the surface effect becomes relatively small and the density profiles 
become almost uniform in each phase. Note again that the electron density is 
different in each case even if we impose the condition ^ = fie in the presence of 
the Coulomb potential. 



In the last of this subsection we summarize the above discussions for the HQ 
phase transition in Table 1. Basically the matter should be treated as a ternary 
system specified by the three kinds of the conserved quantities, baryon number, 
electric charge and lepton number. The proper treatment of the phase transition 
and the mixed phase proceeds by applying the Gibbs conditions. There are some 
kinds of variation instead of the proper treatment: the bulk Gibbs calculation (GCN) 
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Figure 10. Density profiles and the Coulomb potential <p within ID (slab) for ub = 
2.5 no at T = 10 MeV. Here, the neutrino fractions are set to Y Ue = 0.01 (left panel) 
and Y Ue =0.15 (right panel). The cell sizes are Rw = 30 fm in these figures. The 
slab size are R = 10.9 fm. (R/Rw) 3 is fixed to be the optimal value at Y Ve = 0.01. 
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Figure 11. Left panel shows the pressure for PNS matter of T = 30MeV and 
Yl = 0.4. The nuclear density is denoted by no- Right panel shows the density 
profile and the Coulomb potential 3uq assuming the rod structure for the same PNS 
matter. 



may be performed by neglecting the finite-size effects. The bulk Gibbs calculation 
has been also done with the use of the assumption of local charge neutrality (LCN). 
One may regard the matter as a binary system since the chemical equilibrium is 
required only for baryon and lepton numbers in this treatment. In the realistic 
treatment of the mixed phase, the finite-size effects are indispensable. When the 
finite-size effects are taken into account, we see the emergence of the inhomogeneous 
structure called pasta. We have seen that the proper treatment based on the Gibbs 
conditions gives almost the same EOS as the one given by the bulk Gibbs calculation 
with LCN in PNS matter. This is because the scale of inhomogeneity becomes large, 
e.g. Rw — > oo, and an amorphous state is formed, where the finite-size effects are 
irrelevant and local charge neutrality is also attained. This is due to the mechanical 
instability by the neutrino trapping effect. Thus we can say the ternary system is 
effectively reduced to a "binary" system. We attach "binary" in the table to stress 
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this situation. 

The situation somewhat changes for various cases. Neutron star (NS) matter is 
a cold crystallized matter and includes no neutrino. So it should be treated as a 
binary system specified by the baryon number and electric charge. We have seen 
that it is effectively reduced to a "pure" system due to the mechanical instability for 
sufficiently large surface tension or high temperature; the scale of the inhomogeneity 
becomes large again and local charge neutrality is approximately achieved. Conse- 
quently the EOS resembles the one given by the Maxwell construction. The NS-NS 
mergers corresponds to the latter case, and we named "merger (MG)" matter. In 
this sense we listed the system as "pure" in the table for these cases. Note that 
there is one important difference between these cases and PNS matter, while the 
situation looks similar to each other; there is no amorphous state in these cases and 
appears a density jump inside compact stars instead. 

Table 1. Comparison of conditions for the HQ phase transition. 
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3. Possibility to Probe the Properties of the Mixed 
Phase 

As mentioned above, possible existence of the HQ mixed phase in high density region 
can be theoretically expected. However, it might be quite difficult on the Earth 
to reveal the properties of such an extreme state. The one of the possibilities to 
overcome this difficulty may be the direct observation of neutron stars. In fact, since 
the stellar structure of neutron stars strongly depends on their interior properties, 
one could obtain the information of the stellar interior via the observations of stellar 
mass M, radius R, and compactness defined as M/R. 

Once the EOS of matter is obtained, we can extract the density profile (den- 
sity as a function of radial distance from the center) and consequently the mass 
and the radius of compact stars. This is done by numerically solving the Tolman- 
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Oppenheimer-Volkoff (TOV) equation [49] 



dp 
dr 
dm 
dr 



Gme (1 + p/e) (l + 47rr 3 p/m) 



(25) 



r 2 1 — 2Gm/r 

4irr 2 e , 



(26) 



being G the gravitational constant. Starting with a central density e(r = 0) = e c , 
one integrates out Eqs. ([25]) and ([26]) until the surface pressure becomes to zero. 
This gives the stellar radius R and its gravitational mass M = m(R). 

We show the effect of the HQ phase transition on the relation between the EOS 
and bulk structures of compact stars in Fig. [12] The left panel is the EOSs and the 
right panel shows the mass-radius relation of calculated neutron stars. As for the 
EOS, we shows three types of the HQ phase transitions; e.g. the EOS by the bulk 
Gibbs calculation, the Maxwell construction, and full calculation considering the 
finite-size effects. Here, the value of the surface tension is set as a = 40 MeV-fm -2 . 
We also show the uniform hadron and quark matter for comparison. Using them, 
we can calculate the mass-radius relations shown in the right panel. It is well-known 
that the typical masses of neutron stars are distributed in a narrow band close to 
1.4 times the solar mass Mq [50J. The maximum mass of neutron star is enough 
larger than these masses. With inclusion of hyperon degree of freedom, however, 
the softening of matter reduces the possible mass of neutron stars as low as 1.3 Mq, 
which contradicts with the most observations. The blue curves show the cases with 
the HQ mixed phase. The upper most one and the bottom are the results with 
the EOS obtained by the Maxwell construction and the bulk Gibbs calculation, 
respectively. The result with inclusion of the full pasta structures with a = 40 
MeV ■ fm~ 2 lies just below the Maxwell construction case. 

In any case with the HQ mixed phase, the calculated maximum mass is slightly 
above 1.4M0, which might be tolerable in confrontation with the observation. What 
is important is that the mixture of hyperons in nuclear matter softens the EOS, 
which causes the contradiction with the observation, i.e. the "maximum-mass prob- 
lem". By considering the deconfinement transition to quark matter, the EOS gets 
stiffer at higher density region. This stiffness at high density sustains the neutron 
stars with M ~ 1.4M . 

Recently Demorest et.al. have reported that there is a two-solar-mass neutron 
star through the measurement of the Shapiro delay [51]. Motivated by this obser- 
vation, some authors studied the possibility of hard EOS considering exotic matter; 
e.g. Ref.[52j for hadron matter with hyperons, and Ref . |37j the matter including 
quarks. We are also planning to update our EOS models to other models to satisfy 
the observation. 

Another possibility is the observations of gravitational waves radiated from neu- 
tron stars or the observations of stellar oscillations itself, where the observed spec- 
tra can reveal the interior properties |53[ 154] . This technique is well-known as 
"gravitational-wave asteroseismology" or just "asteroseismology," which is similar 
to the seismology in the Earth and the helioseismology in the Sun. In this sec- 
tion, we focus on the possibility to see the properties of the HQ mixed phase via 
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Figure 12. Mass-radius relation of calculated neutron stars. Green curve indicates 
the result with only nucleons.Red curve is that with hyperons. Blue curves are 
results with hadron matter and quark matter by different treatments of the mixed 
phase Note that the Maxwell construction is incorrect in the present case, but it is 
still useful as an eye guideline. 



the observations of gravitational waves in subsection 13.1.1 and stellar oscillations in 
subsection 13. 2. i 

In the last part of this section, we introduce another possibility to know the 
effects of the HQ mixed phase on compact-star phenomena; thermal cooling of 
neutron stars. A cooling theory of neutron stars has been well-discussed but a still 
uncertain issue |55} 156] . Cooling depends on the internal state which is in the range 
of ultrahigh density around/above the normal density. It may be directly connected 
to the EOS discussed in the previous section. 

3.1. Gravitational Waves 

The gravitational waves are oscillations of spacetime itself. Due to their strong 
permeability, the gravitational waves can carry away the raw information about the 
wave sources. So, observing the gravitational waves emitted from neutron stars, it 
could be possible to collect the astronomical data [53], [54l [57] , to reveal the properties 
of dense matter [581 E3 EH EH E2] E3], and to probe the gravitational theory in 
the strong gravitational field [Ml EH ESI E3 EH]- Now, in order to directly detect 
the gravitational waves, several ground-based detectors are in operation, and the 
projects to build the next generation detectors are also moving ahead [69] [70] [71], 
which could enable us to observe the gravitational waves in the not-so-distant future. 

Considering the gravitational waves radiated from the spherically symmetric 
stars, their oscillations can be classified into two families with their parities, i.e., 
odd (axial) and even (polar) parities. Since the odd parity oscillations are com- 
pletely decoupled with the even parity oscillations, one can examine each type of 
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oscillations separately. The odd parity oscillations are incompressible motion, while 
the even parity oscillations involve the deviation of density. Thus, from the ob- 
servational point of view, the even parity oscillations might be important in the 
emission of gravitational waves. Here, in order to determine the specific oscillation 
frequencies in neutron stars, we adopt the relativistic Cowling approximation, where 
the perturbation of gravitational potential should be neglected. Namely, the fluid 
would oscillate on a fixed background spacetime. In practice, even with this simple 
approximation, one can qualitatively see the features for oscillation frequencies of 
emitted gravitational waves. 

During the stellar oscillations, many kinds of gravitational waves could be ra- 
diated. The fundamental (/) and the pressure (p) modes can be excited in the 
spherically symmetric stars without the density discontinuity inside the star, while 
the gravity (g) modes can be also excited with the density discontinuity. As shown 
in Fig. [T3l the stellar models with the EOSs adopted in this article except for the 
EOS with the Maxwell construction have only / and p modes, where the details how 
to determine the specific oscillation frequencies can be seen in Ref. [63J. On the 
other hand, in the stellar models constructed with the EOS adopted the Maxwell 
construction whose masses are more than O.7M , the g mode is also excited as 
well as / and p modes. It should be noticed that since the stellar models with 
M < 0.7Mq do not have density discontinuity inside the star because of the small 
central density, the g mode can not be excited. Since, as mentioned before, the g 
mode can be excited due to the existence of density discontinuity inside the star, 
the observations of the g mode can tell us such an information. In fact, the typical 
frequencies of g modes are in the range from a few hundreds Hz up to kHz, it could 
be possible to detect via the ground-based gravitational wave detectors. Addition- 
ally, as shown in Fig. \13\ the frequencies of g mode for the stellar models adopted 
in this article are almost independent of the stellar mass, i.e., ~ 1.73 Hz. While, as 
shown in Fig. [Til it is found that the g mode frequencies can be well expressed as a 
function of stellar compactness, such as ujM = 0.3130(M/i?) +0.0103. That is, via 
the detailed observations of g mode gravitational waves, one could know the stellar 
properties with the help of the additional observations about the stellar parameters 
such as the mass. 

Furthermore, it is well-known that the / mode frequency can be associated with 
the stellar average density defined as (M/i? 3 ) 1 / 2 [53l[54]. This could be physically 
explained by considering the relation between the sound speed and the propagation 
time of the fluid perturbation inside the star. The / mode frequencies calculated for 
the stellar models with the adopted EOSs are shown in Fig. [15j Roughly speaking, 
one can observe from this figure that the frequencies for the stellar models without 
the density discontinuity, i.e., except for the EOS with the Maxwell construction, 
are almost similar behaviors. Thus, it might be difficult to reveal the properties of 
the HQ mixed phase via the observations of / mode gravitational waves. However, 
it should be emphasized that the / mode frequencies for the stellar model with 
the density discontinuity (with EOS adopted the Maxwell construction) can behave 
qualitatively different from the others. That is, via the observations of not only 
g modes but also / modes gravitational waves, one might be able to prove the 
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Figure 13. The frequencies of / and p\ modes for the stellar models with (a) nucleon 
and (b) the HQ mixed phase adopted the Maxwell construction as a function of the 
stellar mass. In the right panel, the additional frequencies of g modes are also shown 
as well as the frequencies of / and p\ [63] . 
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Figure 14. The normalized eigenfrequencies of g modes with EOSs considering the 
Maxwell construction as a function of the stellar compactness M/R [63J. 

existence of the density discontinuity inside the star. Furthermore, from Fig. \W\ 
for the normalized eigenfrequencies of / and p\ modes, one can obviously see the 
difference between the frequencies of / mode gravitational waves with and without 
density discontinuity as mentioned the above, while the p\ mode gravitational waves 
seem to depend on the adopted EOSs even if the density discontinuity does not exist 
inside the star. Therefore, the detailed observations of p mode gravitational waves 
could reveal the properties of the HQ mixed phase. 

3.2. Shear Oscillations 

Unlike the gravitational waves emitted from neutron stars, there exist the observed 
evidences of the stellar oscillations, i.e., the quasi-periodic oscillations (QPOs) in 
giant flares radiated from soft-gamma repeaters (SGRs). Up to now, three giant 
flares have been observed, such as SGR 0526-66 in 1979, SGR 1900+14 in 1998, 
and SGR 1806-20 in 2004, and the existence of a lot of QPO frequencies in the 
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Figure 15. The frequencies of / modes for the stellar models with different EOSs 
as a function of the stellar average density (M/i? 3 ) 1 / 2 . The results with nucleon, 
hyperon, HQ mixed phase adopted the bulk Gibbs conditions (GCN) or the Maxwell 
construction are shown in the lift panel, while the right panel focuses on the results 
with the EOS considered the finite-size effects. Here the units of each a are MeV • 
fm~ 2 O. 



afterglow of giant flares are found via timing analysis, whose frequencies are in the 
range of tens Hz up to kHz [72]. In order to theoretically explain the observed QPO 
frequencies, many examinations have been done in terms of the torsional oscillations 
in the crustal region and/or the magnetic oscillations (e.g., [751 EH E3 [TBI [TTJ [751 
P79| 1501 EE]). If the observed QPO frequencies are associated with the torsional 
oscillations in neutron star crust, one might be able to obtain the information about 
the crust properties [521 E31 EH ESI 86J. The torsional oscillations are incompressible, 
i.e., odd parity oscillations, where the elasticity plays a role as a restoring force, and 
such elasticity can be characterized by the shear modulus \i [87]. Although the 
derivation of [i for the complex structure have never been done, the expression of ji 
due to the Coulomb interaction in the body center cubic (bcc) in neutron star crust 
is known as 

H = 0.1194ni(Ze) 2 /a, (27) 

where rii, Ze, and a are ion number density, charge of ion, and average ion spacing 
defined as 47ra 3 /3 = 1/rij, respectively |88j . Note that, averaging the all directions 
and assuming the zero temperature, this expression is derived from the Monte Carlo 
calculations. 

On the other hand, as mentioned the above, the possibility to appear the nonuni- 
form structure (the HQ pasta phase) in the high density region inside the neutron 
stars is suggested [TU [19] , which is similar structure to the pasta phase in the neu- 
tron star crust [HE]. If this suggestion is true, similar to the neutron star crust, one 
might be able to consider the shear in the HQ mixed phase, and the torsional oscil- 
lations could be excited even in the HQ mixed phase [89, 90J. In order to examine 
such torsional oscillations, for simplicity, we consider the shear only for the quark 
spherical droplet in the uniform hadronic matter, although the different nonuniform 
structures could appear in the HQ pasta phase, such as cylindrical and slab struc- 
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Figure 16. The normalized eigenfrequencies of / modes (left panel) and pi mode 
(right panel) as a function of the stellar masses. The marks in the figures are the 
same meaning as in Fig. [15] |63j. 

tures. Due to this assumption, the shear in the HQ mixed phase effectively makes 
an underestimation. So, the calculated frequencies of torsional oscillations could be 
corresponding to the lower limit, because the frequencies of torsional oscillations 
can be proportional to the shear speed defined as v s = (fi/p) 1 ^ 2 . Moreover, to adopt 
the expression for fj,, Eq. (|27|) . we consider that rij should be the number density 
of quark spherical droplet in the hadronic sea, while Ze should be the total charge 
included in the quark spherical droplet. As a typical value of the surface tension 
a, we adopt three values, such as a = 10, 20, and 40 MeV • fm" 2 . The detail for 
calculations can be seen in |90j . 

The calculated frequencies of the 1 = 2 fundamental torsional oscillations in the 
HQ mixed phase, o^2) axe shown in Fig. [T7]as a function of the stellar mass. Since 
the frequencies of fundamental torsional oscillations in the neutron star crust are 
around tens Hz [T6|, 187] . one can observe around 10 times larger frequencies in the 
HQ mixed phase. This is why the shear modulus in the mixed phase calculated 
from Eq. (|27p becomes about 10 3 times larger than that in crust. As a result, the 
propagation time with the shear speed becomes about 10 times smaller than that 
in the crust. Then, one can estimate that the frequencies of torsional oscillations in 
the mixed phase could become roughly 10 times as large as those in crust region. 
Additionally, one can see from Fig. [T7] that the frequencies of fundamental torsional 
oscillations depend strongly on the value of a. In practice, the frequencies for a = 20 
and 40 MeV- fm~ 2 are ~40% and ~120 % larger than those for a = 10 MeV- fm~ 2 . 
Thus, if one would identify the observed frequencies as the torsional oscillations 
in the HQ mixed phase, one can probe the properties of such exotic structure. 
Additionally, with the fixed stellar mass, we plot the frequencies of fundamental 
torsional oscillation with £ = 2 in the HQ mixed phase as a function of surface tension 
a in Fig. [THJ From this figure, one can obviously observe the linear relation between 
the frequencies and surface tension. Thus, with the help of the other observation 
of stellar mass, one can make a constraint in the value of a via the observation of 
the frequencies of fundamental torsional oscillations in the HQ mixed phase. Now, 
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Figure 17. With a = 10, 20, and 40 MeV- fin" 2 , the frequencies of the £ = 2 
fundamental torsional oscillations o^2 ar e shown as a function of stellar mass M/Mq 
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Figure 18. Frequencies of fundamental torsional oscillations with I = 2 as a function 
of surface tension with different values of stellar mass |90| . 

we should remark the possibility of the detection of this type oscillation. In fact, if 
the star is exactly spherically symmetric, one might have no chance to observe this 
oscillation, because the HQ mixed phase is surrounded by the hadronic fluid matter 
and the oscillation is completely confined in the HQ mixed phase. But, considering 
the realistic stellar model, there exist the rotational and/or magnetic effects, which 
make a deformation of stellar shape. Under this situation, the torsional oscillations 
with axial parity can be coupled with the even parity oscillations mentioned in Sec. 
13.1.1 Probably, this is one opportunity to observe the torsional oscillation in the 
HQ mixed phase. Anyway, the examinations in the more complicated system will 
be done somewhere. Furthermore, since the resulting frequencies of fundamental 
torsional oscillations are order of 100 Hz, some of the QPO frequencies observed in 
giant flares, for example 150 Hz and even 626.5 Hz in SGR 1806-20 and 155 Hz in 
SGR 1900+14 |72j . might be associated with the torsional oscillations in the HQ 
mixed phase. 



26 



Nobutoshi Yasutake,Tsuneo Noda et al. 



3.3. Cooling of Neutron Stars with the Hadron-quark Mixed Phase 

Now we focus on the cooling theory of isolated neutron stars. Isolated neutron stars 
do not have internal heat source (someone says magnet-breaking heating affects to 
the cooling, but the effect is not so large), and they are born at the sites of super- 
novae explosion where the temperature rises to 30-50 MeV in the core of PNS as 
we described in the last section. Isolated neutron stars evolve thermally by emit- 
ting energy and losing their temperature immediately, as a result the temperature 
becomes to a few MeV (or ~ 10 10 K) in a minute [9T| I92j . In this subsection we 
discuss the cooling of neutron stars after a few minutes from their birth. Hence 
initial temperature is assumed as ~ 1 MeV in our calculations. 

The most efficient energy-release process in neutron stars is the neutrino emis- 
sion, based on neutron beta-decay. The well-known process named direct URCA 
process is combination of neutron beta-decay and inverse reaction. 

n — > p + e~ + u e , p + e~ — > n + v e (28) 

Working forward and inverse (electron capture) process, electron neutrinos and anti- 
electron neutrinos are emitted and carry some energy. Emitted neutrinos travel 
through the stars without any interactions due to the large mean free path. However, 
due to the abundances of proton and neutron, direct URCA process is kinematically 
forbidden in the neutron stars [19]. To solve this issue, the modified URCA process 
was promoted, which involves another neutron to conserve momentum. 

n + n — > n + p + e~ + u e , n + p + e~^rn + n + v e (29) 

The modified URCA process has a lower emissivity; it is about five orders of mag- 
nitude weaker than one with the direct URCA. 

Other processes also work on neutron stars cooling, such as Bremsstrahlung 
(neutrino emission) and photon emission, but they work dominantly in the crust 
(Bremsstrahlung), or the surface (photon emission) of neutron stars, and their emis- 
sivities are weaker than the core processes. 

Modified URCA and other processes could explain most of neutron stars obser- 
vational data, with some exception. Therefore the "Exotic (Non-Standard) Cooling" 
which works in some kinds of exotic phases, is introduced and the ordinary cool- 
ing is called the "Standard Cooling" [93J. The standard cooling works in all the 
neutron stars, and the exotic cooling works in some of neutron stars. The exotic 
cooling works above the threshold density, due to some kind of exotic situations 
where the exotic phase appears. The exotic situation can be said as pion condensa- 
tion, kaon condensation, existence of hyperons or quark matter |49j . In the exotic 
state, the neutrino emission process plays as the direct URCA, due to the degrees 
of freedom. Once the exotic cooling process works, star cools faster than standard 
one (eg. [91 [95]). 

Stars with the exotic cooling have larger central density than ones with the 
standard cooling. The mass of stable neutron stars increases with the central density, 
therefore stars with the exotic cooling have larger mass than standard ones [96] . 
Also, the neutrino emissivities of these processes including the modified URCA have 
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density dependent; it is monotonically increasing with the density. It is believed that 
lighter neutron stars cools slower and heavier one cools faster. 

The matter in neutron stars is well degenerated; neutrons and protons can create 
the Cooper pairs and condensate to superfluid (superconducting for protons) state 
at some temperature. The transition between normal to superfluid state, the latent 
heat escapes by causing neutrino pair creation [97]. This process is called as PBF 
(Pair Breaking and Formation) process, and it works only where the temperature 
of the region is near to the critical temperature of superfluidity. In neutron stars, 
there are three possible nucleon superfluid states ( x 5o and 3 P2 for neutrons, and 
1 Sq for protons), and they might have three PBF working regions. The neutron 3 P2 
state works at the highest density of these three, and it may become the dominant 
cooling process at specific density and time [98l [99]. The transition does not last 
long time and occurs just once in a density region, therefore the integrated effect is 
limited. However, once matter becomes to superfluid state, thermodynamical prop- 
erties decreases and the emissivity of neutrino also decreased. Nucleon superfluidity 
works both of acceleration and breaking of cooling. The neutrino emission by the 
transition to the superfluid state works all the neutron stars, and called as "minimal 
cooling" including the modified URCA, Bremsstrahlung and PBF [100 . 

There are some cooled stars whose effective temperature cannot be explained 
by the standard cooling scenario, which does not include the neutrino emission 
caused by nucleon superfluidity. It needs stronger cooling process as in the case 
of J0205+6449 in 3C58 (hereafter, "3C58") or Vela pulsar (B0833-45). Also an 
accreting neutron star SAX J1808 requires strong cooling. 3C58 and Vela may be 
explained by the minimal cooling model which includes nucleon superfluidity |101j . 
However, SAX J1808 needs stronger cooling than the minimal cooling [102] . 

In recent observation, two important data for cooling of neutron stars are ob- 
served and analyzed. One is the central source of Cassiopeia A (hereafter "Cas 
A") [10311104] and the other is J1614-2230 [5l]. They give strong constraints on the 
cooling of neutron stars. J1614-2230 is discussed in the beginning of this chapter, 
we introduce about the Cas A issue here. 

Cas A is the youngest-known supernova remnant in our Milky Way and it locates 
~ 3.4 kpc distant from the solar system [105] , The supernova explosion occurred 
about 330 years ago, but due to absorption by interstellar medium, there are no 
exact historical records except for unclear detection by J. Flamsteed in 1680 [106 . 
Recently, |103] and [104] have analyzed the X-ray spectra of Cas A. They give the 
effective temperature and possible regions occupied by the mass-radius relations. 
Since Cas A is the isolated remnant, uncertainty of the mass-radius relation could 
be large. The lowest mass obtained from x 2 fitting is about 1.5M . Considering the 
age of t = 330 yr, the effective temperature T e s of Cas A must occupy a point of a 
cooling curve due to the standard cooling scenario on the (T e g — t) plane. This gives 
strong constraint on the EOS and cooling processes. Furthermore, [104J reported 
the observation of T e g for Cas A in recent 10 years. [107] . [108] and [109] insist 
that the rapid decrease in T e fj during the years shows that the transition to nucleon 
superfluidity occurs. 
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3.3.1. Cooling of Hybrid Stars with the Hadron-quark Mixed Phase 

Considering the HQ mixed phase for the EOS of neutron star, it is important 
to apply for the cooling theory. The quark matter is known as an exotic matter 
which causes strong neutrino emission. Therefore, once the quark phase appears, 
strong neutrino emission cools the star quickly. The neutrino emissivity is calculated 
as [TTO] 

e vq = 8.8 x I0 26 a s — y e 1/3 T 9 6 erg s" 1 cm- 3 . (30) 
n 

This emissivity can work inside of quark matter, and we show the following assump- 
tion. 

We constructed a model which includes the HQ mixed phase. Considering the 
first order phase transition between the hadron and quark phases, it would be plau- 
sible that both phases coexist and form some kinds of the mixed phase. It has been 
shown that the mixed phase could form geometrical structures [18]; [19] have made 
EOS of mixed-phase under a Wigner-Seitz (hereafter "WS") approximation using 
the MIT Bag model for the quark phase at finite temperature. 

We employ an EOS with the same framework using the bag constant B = 
100 MeV fm -3 , the coupling constant as = 0.2, and the surface tension param- 
eter a = 40 MeV fm -2 . For the hadron phase, we adopt results of the BHF theory 
including hyperons, A, and E~ [23j [231 113 [22] • Since the BHF results is inappro- 
priate for low-density matter in the crust, we apply EOS of BPS [lllj for the crust. 
The EOS gives the maximum mass 1.53Mq with the radius 8.6 km, and the mass 
lies within the limits of the observation of Cas A. 

Using the WS approximation, we obtain a cell radius of each phase, and calculate 
the volume fraction of quark matter in the mixed phase as seen in Fig. [19j 

It is difficult to calculate the neutrino emissivity in the mixed phase directly. 
Therefore, the volume fraction F is multiplied to the original quark neutrino emis- 
sivity e Vtq of (i30j) : The total emissivity by quarks is set to be e v = Fe V) g. We adopt 
the well-known standard cooling process for hadronic matter [93] : Modified URCA 
process for the higher density region and Bremsstrahlung process for the crust. 

Considering the HQ mixed phase, it becomes possible to reduce the neutrino 
emissivity by quarks. However, the cooling rate is still too strong to explain cooled 
stars with lower limit (such as Vela) , it requires some kind of reduction of neutrino 
emissivity. 

3.3.2. Cooling of Hybrid Stars with the Hadron-quark Mixed Phase and 
Color Superconductivity 

Fortunately, there is possible factor of reduction of neutrino emissivity in the quark 
phase. It is the color superconducting (hereafter "CSC") phase of quarks |112j . 
There are some sorts of quark pairings such as CFL or 2SC according to the degrees 
of freedom of quark flavor and color. In the CFL phase, all colors and flavors quarks 
have a common Fermi momenta, and all quarks create pairs and condensated. Mean- 
while u and d quarks of two colors have same Fermi momenta in the 2SC phase, 
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Figure 19. Volume fractions of quark matter phase having particular geometrical 
structures with a bag constant B = 100 MeV fm -3 and a coupling constant a s = 

0.2 pans]. 




log f [yr] 



Figure 20. Cooling curves with the HQ mixed phase. Solid and dashed curves 
denote cooling curves of the hybrid stars with the HQ mixed phase, and the dotted 
curve denotes the cooling curve of the neutron star with standard cooling. The 
curves show the models of hybrid stars with the HQ mixed phase become too cold 
to explain observed data. 
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and some quarks are unable to create pairs [3]. The CFL phase is likely to appear 
at higher density than the 2SC phase |113j . It is considered that the energy gap 
of A ~ 10 MeV is very large compared with the temperature of the center of com- 
pact stars , Tq ~ keV,(eg. |114| ). Once matter becomes superconducting, neutrino 
emissivity must be suppressed similar to the hadronic superfluidity |115j . Due to 
the large energy gap and it could be proportional to exp (— A/feT), where T is the 
temperature at the relevant layer and /cb is the Boltzmann constant. Therefore, in 
the color superconducting phase with a large energy gap (A > T), neutrino emis- 
sivity by quarks is almost negligible [3J. We note that we do not need to consider 
which kinds of CSC pairing appears, if we adopt a large enough energy gap. 

We assume that the energy gap in the color superconducting phase is very large 
and the phase appears in the density region above the threshold density px- For 
simplicity, we use the critical volume fraction of the quark matter in the HQ mixed 
phase Fq instead of px- If the matter has a volume fraction F > Fq at the density 
in the layer of the mixed phase, quarks become the color superconducting state. 

We select stellar masses of compact stars to be 1.0, 1.3, and 1.5Mq whose central 
densities are 1.48, 1.82, and 2.67 xlO 15 g cm -3 , respectively. We take the critical 
fraction Fq to be 0.1, 0.125, and 0.2 which are appropriate to explain the observa- 
tions with our scenario. Results are shown in Fig. [2X1 with available observational 
values. Most of the observational data are taken from the Table 7.2 in |116j : |117j 
for 3C58, P3H] for Crab, pTS] for Vela, [120] for 0656+14, [m] for 1706-44, pT6] 
for 1811-1925, [T22] for 1823-13, [US] for 2334+61. Other data are taken from [124] . 
The data of Cas A is attached |104j as the most youngest compact star. The cool- 
ing curves transit from hotter region to cooler region for the parameter between 
0.1 < Fq < 0.2. As indicated in the middle panel of Fig. [2T| cooling curves split 
into two regions for Fq = 0.125 and stars with larger masses cool slower than those 
with lighter masses. Since lighter mass stars cool faster, they suit for the 3C58 case 
who does not have lower limit of the effective temperature observation. However, 
calculated cooling curves are inconsistent with the Vela case which has lower limit 
and the quark cooling is still too strong to explain this case. 

Since the neutrino emissivity of the quark phase involves large uncertainty, we 
have calculated the additional cooling curves for the mass 1.03M Q in case of Fq = 
0.125 with the neutrino emissivity reduced by a factor of 0.1 and 0.01. There are 
some possible factors of this reduction for neutrino emissivity accompanying quark 
/3-decay, such as increase of the abundance of strange quarks; decreasing in electron 
numbers inside MP, leads to the reduction of the neutrino emissivity |110j . The 
presence of 2SC at low density also suppress the emissivity; [18] discussed that 
the abundance of quarks in the mixed phase changes and it may cause the color 
superconducting phase. We suppose that the reduction of the emissivity originates 
from the above physical processes. If the emissivity of quarks is reduced by these 
factors, the observation of Vela can be explained as shown in the middle panel of 
Fig. EU Also, we calculated the effect of PBF and the carbon envelope with our 
model to fit the recent Cas A data [104] . and the result is shown in the thick gray 
line in the middle panel of Fig. [2TJ 

Considering the HQ mixed phase, the cooling theory of neutron star requires the 
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Figure 21. Cooling curves with color superconducting quark phases |125j . Solid, 
dotted, and dashed lines denote the models with the masses of 1.50, 1.32, and 
1.03 Mq, respectively. Thick gray line on the middle panel denotes a 1.5OM0 model 
with PBF and carbon envelope. Dot-dashed lines with marks on the middle panel 
indicate the model of the mass 1.O3M0 except for the neutrino emissivity in normal 
quark phase multiplied by 1/10 and 1/100 for the lines with marks of circle and 
triangle, respectively. 



existence of some kind of color superconducting phase with large gap. If the mass 
of an isolated neutron star other than Cas A is detected in near future, it might 
help to probe the existence of color superconducting phase, quark matter and the 
HQ mixed phase. 



4. Summary and Concluding Remarks 

We have discussed features of the HQ mixed phase in compact stars, where multi 
components are in chemical equilibrium. Such phase transitions are generally not 
congruent and the Maxwell construction, which is a common method to obtain the 
equation of state (EOS) of matter with a single component, cannot be applied any 
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more and one has to solve the Gibbs conditions instead. Using the bulk calculation, 
i.e., without considering the surface tension and the Coulomb interaction among two 
phases, we have presented basic concepts of the non-congruent transition. Although 
these concepts are helpful in discussing and analyzing the properties of the phase 
transitions in compact stars, we must bear in mind that they are obtained in an 
idealized and simplified situation. For realistic cases, we have to take into account 
rather complicated finite-size effects such as the surface tension and the Coulomb 
repulsion. 

When the finite-size effects are present, crystalline structures of regular geome- 
tries called "pasta structures" emerge in the mixed phase. Charge screening, i.e., 
the rearrangement of charged particle densities in the presence of the Coulomb in- 
teraction is also emphasized; it sometimes causes an instability of the geometrical 
structures in the mixed phase. In the extreme case the EOS resembles the one 
given by the Maxwell construction, and uniform nuclear matter can exist in the 
mechanically unstable region. 

We note that the maximum mass of hybrid stars in our calculation is rather low 
in the light of a recent observation [51j . since contamination of hyperons considerably 
softens EOS. We have found that the transition of hadronic matter to quark matter 
suppresses the appearance of hyperons, but the maximum mass is still low due to 
the rather soft EOS of quark matter. Therefore, we need other EOSs of hadron and 
quark matter to circumvent this situation. 

So far, most works on the pasta structures have used the WS cell with ansatz 
about the geometrical structures like droplet, rod, slab and so on. If we want to 
see the possibility of the intermediate shapes, or study the mixed phase in a more 
realistic way, Frameworks without the WS approximation or ansatz on the geo- 
metrical structures are desirable. Recently we have performed a three-dimensional 
calculation of non-uniform nuclear matter based on the relativistic mean-field model 
and the Thomas- Fermi approximation |12j . Introducing a cubic cell with periodic 
boundary conditions and discretizing it into grids, we numerically solve the coupled 
field equations for meson mean-fields and the Coulomb potential. Randomly dis- 
tributing fermions (n,p,e) on the grid points as an initial condition, we relax their 
density distributions to attain the uniformity of the chemical potentials. We have 
shown that typical pasta structures, obtained in the previous studies with the WS 
approximation, also appear in the three-dimensional calculation [12J. 

Possibilities of observing the signals of the mixed phase has been suggested in 
the spectra of the gravitational waves [63] • We have found that one could know 
the existence of density discontinuity inside the star via observing the gravitational 
waves of not only g mode but also / mode. The detailed observations of p mode 
gravitational waves could reveal the properties of the HQ mixed phase because 
the p\ mode gravitational waves depend on the adopted EOSs even if the density 
discontinuity does not exist. 

QPO in giant flares from SGR is another probe from neutron stars. If the 
observed QPO frequencies are associated with the torsional oscillations in neutron 
star crust, one might be able to obtain the information about the crust properties. 
In fact, we have pointed out that some of the QPO frequencies observed in giant 
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flares might be associated with the torsional oscillations in the HQ mixed phase. 

Phenomenological implications of pasta structures should deserve more elaborate 
studies of their thermodynamical or dynamical properties, e.g., the neutrino opacity 
or the thermal conductivity may be affected by the pasta structures, and their 
elasticity may cause a discontinuous change of internal structure of compact stars. 
In particular, the crystalline structure of the HQ mixed phase in the core region 
may bring about novel dynamical effects on compact-star phenomena by way of 
deformation or oscillation. 

The other important application of the EOS with the HQ mixed phase is the 
cooling theory of neutron stars, the long time issue discussed for decades. The 
cooling of neutron stars strongly depends on the neutrino emissivity of the core of 
stars. Simply applying the HQ mixed phase with normal quark matter, stars cool 
rapidly and do not match observations. However, including color superconductivity, 
the neutrino emissivity of the quark phase is much suppressed, and the cooling of 
stars with the HQ mixed phase does not conflict with observations. 

Taking account of a color superconductivity in the HQ mixed phase, it became 
possible to have heavier stars cooling slower and lighter ones faster. We cannot get 
this situation without considering the HQ mixed phase. If one observes the masses 
of isolated neutron stars and the mass-temperature distribution can be described 
by the above scenario, it proves the existence of the color superconductivity in the 
HQ mixed phase. 
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